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The  Relationship  Between 
Frictional  Resistance  and 
Roughness  for  Surfaces 
Smoothed  by  Sanding 

An  experimental  investigation  has  been  carried  out  to  document  and  relate  the  frictional 
resistance  and  roughness  texture  of  painted  surfaces  smoothed  by  sanding.  Hydrodynamic 
tests  were  carried  out  in  a  towing  tank  using  a  flat  plate  test  fixture  towed  at  a  Reynolds 
number  (ReL)  range  of  2.8 X 10°  —5.5 X Hr  based  on  the  plate  length  and  freestream 
velocity.  Results  indicate  an  increase  in  frictional  resistance  coefficient  (CF)  of  up  to 
7.3%  for  an  unsanded,  as-sprayed  paint  surface  compared  to  a  sanded,  polished  surface. 
Significant  increases  in  CF  were  also  noted  on  surfaces  sanded  with  sandpaper  as  fine  as 
600-grit  as  compared  to  the  polished  surface.  The  results  show  that,  for  the  present 
surfaces,  the  centerline  average  height  (Ra)  is  sufficient  to  explain  a  large  majority  of 
the  variance  in  the  roughness  function  (A  U  + )  in  this  Reynolds  number  range. 

[DOI:  10.1115/1.1459073] 


Introduction 

Many  practical  engineering  applications  involve  turbulent  flows 
over  surfaces  that  have  been  smoothed  by  sanding.  Examples 
range  from  sailing  hulls  to  models  for  wind  and  water  tunnels. 
While  a  great  deal  of  drag  data  has  been  generated  for  sandgrain 
roughness  (most  notably  Nikuadse’s  experiments  on  monodis- 
perse,  closely-packed  sand  [1]),  there  are  few  reliable  data  for 
sanded  surfaces  in  which  the  surface  is  well  documented.  This  is 
noteworthy  since  sanded  surfaces  form  a  much  larger  presence  in 
engineering  applications  than  sandgrain  roughness.  The  purpose, 
therefore,  of  the  present  investigation  is  to  study  the  frictional 
resistance  of  sanded  surface  roughness. 

A  large  body  of  basic  research  has  focused  on  the  effect  of 
surface  roughness  on  frictional  resistance.  Hama  [2],  Ligrani  and 
Moffat  [3],  Krogstad  and  Antonia  [4]  and  many  others  have  in¬ 
vestigated  the  structure  of  the  turbulent  boundary  layer  over  rough 
surfaces.  Raupach  et  al.  [5]  give  a  review  of  much  of  this  work. 
Studies  focusing  on  the  frictional  resistance  of  ship  bottom  paints 
have  also  been  made.  Grigson  [6],  Townsin  et  al.  [7],  Musker  [8], 
and  Lewkowicz  and  Musker  [9]  have  all  investigated  these  sur¬ 
faces,  and  their  results  indicate  that  as-sprayed  antifouling  coat¬ 
ings  have  significantly  higher  frictional  resistance  than  a  smooth 
surface.  An  entire  workshop  was  devoted  to  the  subject  of  ship 
hull  roughness  and  drag  [10].  However,  most  of  the  research  into 
the  drag  on  marine  paints  has  centered  on  predicting  the  economic 
penalty  of  hull  roughness  on  commercial  ships,  where  sanding  is 
unfeasible. 

The  mean  velocity  profile  in  a  turbulent  boundary  layer  can  be 
expressed  as: 

t/+  =  — ln(v  +  )  +  6  +  2<n(y/£>)n//r  (1) 

K 

Clauser  [11]  noted  that  for  rough  wall  flows,  the  log-law  intercept 
is  shifted  downward  and  that  the  shift  correlates  with  k+ ,  the 
roughness  Reynolds  number,  defined  as  the  ratio  of  the  roughness 
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length  scale,  k,  to  the  viscous  length  scale,  vlU T.  This  downward 
shift,  AU+ ,  called  the  roughness  function,  can  be  used  to  express 
the  mean  velocity  profile  for  rough  wall  flows: 

U+  =  -\n{y  +  )  +  B-\U+  +  2u){yl8)YllK  (2) 

K 

Hama  [2]  showed  that  by  evaluating  Eqs.  (2)  and  (3)  at  y=  S,  the 
roughness  function  could  be  found  by  subtracting  the  rough  wall 
intercept  from  the  smooth  wall  intercept,  B,  at  the  same  value  of 
Re ,5* .  The  roughness  function  therefore  can  be  expressed  as: 


It  should  be  noted  that  Eq.  (3)  is  only  valid  provided  both  IT  and 
the  velocity  defect  profile  are  the  same  for  the  rough  and  smooth 
walls.  The  experimental  evidence  for  this  is  somewhat  contradic¬ 
tory.  Some  research  indicates  that  surface  roughness  increases  TI 
(e.g.  [4])  and  alters  the  velocity  defect  profile  (e.g.  [12]),  while 
other  studies  (e.g.  [13])  indicate  that  these  are  unchanged  by 
roughness.  In  the  present  study,  no  mean  velocity  profiles  were 
made.  Therefore,  the  determination  of  the  roughness  function  re¬ 
quired  the  explicit  assumption  that  II  and  the  velocity  defect  pro¬ 
file  are  unchanged  by  the  roughness  to  be  made.  Future  studies  are 
planned  that  will  include  measurements  of  the  mean  velocity  pro¬ 
files  over  these  sanded  surfaces. 

A  universal  roughness  function  for  a  given  class  of  surfaces  can 
be  defined  if  k  is  related  directly  to  the  surface  profile.  Nikuradse’s 
[1]  pipe  flow  experiments  on  closely-packed,  uniform  sand  rough¬ 
ness  show  that  this  roughness  type  has  a  universal  roughness  func¬ 
tion  with  k  simply  being  the  diameter  of  the  individual  sandgrains. 
The  results  from  Nikuradse’s  experiments  have  been  used  to  ex¬ 
plain  the  behavior  of  generic,  naturally  occurring  surface  rough¬ 
ness.  This  is  evidenced  by  the  widespread  use  of  the  equivalent 
sandgrain  height,  ks .  This  parameter  is  defined  as  the  sandgrain 
height  in  Nikuradse’s  experiment  that  has  the  same  roughness 
function  in  the  fully-rough  regime  as  the  surface  of  interest.  The 
use  of  ks  is  attractive  because  it  is  simple,  but  is  also  problematic 
because  it  is  not  physically  related  to  the  surface  roughness  profile 
for  generic  surfaces  of  engineering  interest.  Most  naturally  occur¬ 
ring  rough  surfaces  do  not  behave  like  Nikuradse  sand  surfaces. 
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Colebrook  [14]  first  demonstrated  this  in  a  study  of  the  irregular 
surface  roughness  in  pipes  resulting  from  the  manufacturing  pro¬ 
cess. 

Nikuradse’s  roughness  function  for  uniform  sand  roughness  has 
led  to  the  critical  roughness  height  concept.  This  concept  states 
there  exists  some  critical  roughness  height  for  surfaces  below 
which  there  is  no  increase  in  drag.  This  is  termed  the  hydrody- 
namically  smooth  condition.  In  this  condition  the  individual 
roughness  elements  are  small  enough  to  be  completely  submerged 
in  the  viscous  sublayer  region  of  the  boundary  layer.  In  order  to 
have  a  hydrodynamically  smooth  surface,  k+  must  be  less  than  a 
critical  value  ranging  from  2.25  to  5.  For  this  surface  type,  if  the 
viscous  length  scale  is  known,  a  critical  roughness  height  may  be 
specified  for  a  surface  below  which  a  reduction  in  roughness 
height  causes  no  concomitant  reduction  in  drag.  A  recent  paper  by 
Bradshaw  [15]  questions  the  existence  of  a  critical  roughness 
height  on  theoretical  grounds.  He  argues  that  the  roughness  func¬ 
tion  should  go  asymptotically  to  0  in  the  limit  as  k+  goes  to  0. 

Granville  [16]  offers  three  alternative  methods  for  determining 
the  roughness  function  of  a  surface  experimentally  using  pipe 
flow,  towed  flat  plates,  and  rotating  disks,  respectively.  The  pro¬ 
cedure  given  for  towed  flat  plates  was  used  in  the  present  inves¬ 
tigation  to  determine  the  roughness  function.  Further  details  are 
given  in  the  Results  and  Discussion  section  of  the  paper.  It  should 
be  noted  that  A  U+  can  also  be  obtained  directly  by  measuring  the 
mean  velocity  profile  over  a  rough  wall.  Once  A.U+ 
=  A U+  (k+ ,[!])  f°r  a  surface  is  known,  it  can  be  used  in  a  bound¬ 
ary  layer  code  or  a  similarity  law  analysis  to  predict  the  drag  of 
any  body  covered  with  that  roughness. 

A  great  deal  of  effort  has  been  made  to  correlate  the  roughness 
function  for  a  surface  with  its  roughness  statistics.  This  would 
allow  the  drag  change  to  be  predicted  based  on  knowledge  of  the 
surface  profile  alone.  However,  development  of  a  universal  rela¬ 
tionship  to  correlate  the  roughness  function  to  the  surface  rough¬ 
ness  length  scales  has  been  illusive.  Several  researchers  have  at¬ 
tempted  to  correlate  the  roughness  function  with  a  roughness 
height  and  density  parameter  for  relatively  simple  uniform  rough¬ 
ness  [17,18].  Koch  and  Smith  [19]  and  Acharya  et  al.  [12]  both 
looked  at  the  effect  of  machined  roughness  on  frictional  resis¬ 
tance.  Acharya  et  al.  found  that  collapsing  the  roughness  functions 
for  these  surfaces  to  a  universal  curve  using  Ra  alone  was  not 
possible  and  suggested  that  the  deviation  in  the  slope  angles  of  the 
roughness  might  allow  better  correlation.  Both  Townsin  et  al.  [7] 
and  Musker  [8]  have  proposed  correlations  that  include  roughness 
height  as  well  as  texture.  Townsin  proposed  that  a  height  param¬ 
eter,  h,  based  on  the  first  three  even  moments  of  the  profile  power 
spectral  density  reasonably  collapsed  the  roughness  functions  for 
ship  hull  surfaces.  Townsin  and  Dey  [20]  give  the  following  for¬ 
mulation  for  the  roughness  function  of  ship  hull  coatings  based  on 
their  modified  roughness  Reynolds  number: 

AC/+  =  -ln(l  +0.18/i  +  )  (4) 

K 

where 

h=  \Jam0m2 


r tt/ls 

1 

J  2ttILp 


E  y"d  y 


(5) 


Musker  [8]  suggests  an  alternative  roughness  scale  that  incor¬ 
porates  the  skewness  and  kurtosis  of  the  roughness  height  distri¬ 
bution.  Grigson  [6]  asserts  that  the  statistics  of  the  surface  profile 
alone  cannot  be  relied  upon  to  predict  the  roughness  function.  He 
contends  that  only  experimental  testing  of  the  surface  of  interest 
allows  accurate  determination  of  the  roughness  function.  It  is  of 


TOWING  CARRIAGE  DIRECTION 


Fig.  1  Schematic  of  the  flat  plate  test  fixture 


note  that  Grigson’s  results  indicate  that  the  roughness  functions  of 
some  ship  hull  coatings  do  not  behave  as  either  Nikuradse  or 
Colebrook-type  functions  and  may  have  multiple  inflection  points. 

The  goal  of  the  present  experimental  investigation  is  to  docu¬ 
ment  the  frictional  resistance  and  surface  roughness  of  a  range  of 
sanded  paint  surfaces.  An  attempt  to  identify  a  suitable  roughness 
parameter  relating  the  physical  roughness  and  the  roughness  func¬ 
tion  for  this  particular  class  of  surfaces  is  made.  The  results  are 
then  scaled  up  to  a  planar  surface  using  the  similarity  law  proce¬ 
dure  of  Granville  [21]  to  predict  the  effect  of  the  present  rough¬ 
ness  on  the  frictional  resistance  of  a  plate  of  the  order  of  length  of 
typical  sailing  vessels. 

Experimental  Facilities  and  Method 

The  present  experiments  were  conducted  in  the  115  m  long 
towing  tank  facility  at  the  United  States  Naval  Academy  Hydro¬ 
mechanics  Laboratory.  The  width  and  depth  of  the  tank  are  7.9  m 
and  4.9  m,  respectively.  The  towing  carriage  has  a  velocity  range 
of  0-7.6  m/s.  In  the  present  investigation  the  towing  velocity  was 
varied  between  2.0  m/s-3.8  m/s  (ReL=2.8X  106-5.5X  106).  The 
velocity  of  the  towing  carriage  was  measured  and  controlled  using 
an  encoder  on  the  rails  that  produce  4000  pulses/m.  Using  this 
system,  the  precision  uncertainty  in  the  mean  velocity  measure¬ 
ment  was  <0.02%  over  the  entire  velocity  range  tested.  The 
working  fluid  in  the  experiments  was  fresh  water,  and  the  tem¬ 
perature  was  monitored  to  within  ±0.05°C  during  the  course  of 
the  experiments  using  a  thermocouple  with  digital  readout.  Fur¬ 
ther  details  of  the  experimental  facility  are  given  in  [22]. 

Figure  1  shows  a  schematic  of  the  test  fixture  and  plate.  The  flat 
test  plate  was  fabricated  from  304  stainless  steel  sheet  stock  and 
measured  1.52  m  in  length,  0.76  m  in  width,  and  3.2  mm  in 
thickness.  Both  the  leading  and  trailing  edges  were  filleted  to  a 
radius  of  1.6  mm.  No  tripping  device  was  used  to  stimulate  tran¬ 
sition.  The  overall  drag  of  the  plate  was  measured  using  a  Model 
HI-M-2,  modular  variable-reluctance  displacement  force  trans¬ 
ducer  manufactured  by  Hydronautics  Inc.  An  identical  force  trans¬ 
ducer,  rotated  90  deg  to  drag  gage,  was  included  in  the  test  rig  to 
measure  the  side  force  on  the  plate.  The  purpose  of  the  side  force 
gage  was  to  ensure  precise  alignment  of  the  plate.  This  was  ac¬ 
complished  by  repeatedly  towing  the  plate  at  a  constant  velocity 
and  adjusting  the  yaw  angle  of  the  test  fixture  to  minimize  the  side 
force.  Once  this  was  done,  no  further  adjustments  were  made  to 
the  alignment  over  the  course  of  the  experiments.  The  side  force 
was  monitored  throughout  to  confirm  that  the  plate  alignment  did 
not  vary  between  test  surfaces.  Both  of  the  force  transducers  used 
in  the  experiments  had  load  ranges  of  0-89  N.  The  combined  bias 
uncertainty  of  the  gages  is  ±0.25%  of  full  scale.  Data  were  gath¬ 
ered  at  a  sampling  rate  of  100  Hz  and  were  digitized  using  a 


Journal  of  Fluids  Engineering 


JUNE  2002,  Vol.  124  /  493 


16-bit  A/D  converter.  The  length  of  the  towing  tank  dictated  the 
sampling  duration.  This  ranged  from  ~30  s  of  data  per  test  run  at 
the  lowest  Reynolds  number  to  ~11  s  of  at  the  highest  Reynolds 
number.  The  overall  drag  was  first  measured  with  590  mm  of  the 
plate  submerged.  This  was  repeated  with  25  mm  of  the  plate  sub¬ 
merged  in  order  to  find  the  wavemaking  resistance  tare.  The  dif¬ 
ference  between  the  two  was  taken  to  be  the  frictional  resistance 
on  the  two  565  mm  wide  by  1.52  m  long  faces  of  the  plate.  The 
tests  were  repeated  ten  times  for  each  surface  and  velocity.  The 
results  presented  are  the  means  of  these  runs. 

A  single  test  plate  was  used  for  all  the  towing  experiments.  This 
was  done  to  ensure  that  any  differences  in  the  drag  measured  were 
due  to  the  surface  condition  of  the  plate  and  not  small  variations 
in  leading  edge  shape,  plate  flatness,  and  other  factors  that  could 
have  varied  between  multiple  test  plates.  The  plate  was  initially 
painted  with  several  coats  of  marine  polyamide  epoxy  paint 
manufactured  by  International  Paint.  This  surface  condition  was 
termed  the  “unsanded”  condition.  After  hydrodynamic  testing,  the 
plate  was  wet  sanded  with  60-grit  sandpaper.  This  surface  was 
referred  to  as  the  “60-grit  sanded”  condition.  Subsequent  to  hy¬ 
drodynamic  testing  under  the  60-grit  sanded  condition,  the  entire 
process  was  repeated  for  the  “120-grit  sanded,”  “220-grit 
sanded,”  “400-grit  sanded,”  and  “600-grit  sanded”  surface  con¬ 
ditions.  After  hydrodynamic  testing  of  the  600-grit  sanded  surface 
the  plate  was  wet  sanded  up  to  1800  grit  and  polished  with  a 
buffing  wheel  using  Maquire’s  swirl  remover  polishing  com¬ 
pound.  This  surface  is  referred  to  as  the  “polished”  condition.  All 
the  sanding  in  the  present  experiment  was  carried  out  by  hand 
with  the  aid  of  a  sanding  block  using  small  circular  motions. 
Prussian  blue  was  used  to  dye  the  surface  before  sanding  with 
finer  grit  paper  to  ensure  the  entire  surface  was  sanded  and  to 
reveal  the  surface  scratches  left  behind  by  the  previous  grit  so 
they  could  be  removed.  The  surfaces  were  cleaned  with  water  and 
a  soft  cloth  between  surface  treatments  to  remove  grit  and  detritus 
left  behind  by  the  sanding  process. 

The  surface  profiles  of  the  test  plates  were  measured  using  a 
Cyber  Optics  laser  diode  point  range  sensor  (model  #PRS  40) 
laser  profilometer  system  mounted  to  a  Parker  Daedal  (model 
#106012BTEP-D3L2C4M1E1)  two-axis  traverse  with  a  resolu¬ 
tion  of  5  /mi.  The  resolution  of  the  sensor  is  1  /im  with  a  laser 
spot  diameter  of  10  /on.  Data  were  taken  over  a  sampling  length 
of  50  mm  and  were  digitized  at  a  sampling  interval  of  25  /im.  Ten 
linear  profiles  were  taken  on  each  of  the  test  surfaces.  A  single 
three-dimensional  topographic  profile  was  made  on  each  of  the 
surfaces  by  sampling  over  a  square  area  2.5  mm  on  a  side  with  a 
sampling  interval  of  25  /un. 

The  roughness  statistics  were  calculated  using  the  linear  pro¬ 
files  from  each  of  the  surfaces.  All  were  calculated  without  using 
a  long  wavelength  cutoff  (effectively  the  cutoff  was  the  sampling 
length,  50  mm)  and  using  25  mm,  10  mm,  5  mm,  2  mm,  and  1 
mm  long  wavelength  filters.  The  highpass  filtering  was  carried  out 
using  a  Butterworth  digital  filter  and  the  long  wavelength  cutoffs 
were  chosen  to  be  in  the  range  used  by  Musker  [8]  and  Townsin 
and  Dey  [20].  The  purpose  of  the  filtering  was  to  remove  surface 
waviness  which  has  little  effect  on  the  drag.  Musker  says  that  the 
long  wavelength  cutoff  should  be  set  equal  to  the  size  of  the  large 
energy-containing  eddies  near  the  surface,  and  he  suggests  using 
the  Taylor  macro-scale.  In  the  present  investigation  no  spatial  tur¬ 
bulence  correlations  were  available  from  which  to  calculate  the 
Taylor  macro-scale,  so  roughness  statistics  using  a  range  of  long 
wavelength  cutoffs  were  calculated. 

Some  of  the  roughness  statistics  calculated  for  the  surfaces  in¬ 
cluded  the  centerline  average  height,  Ra  ,  given  as: 

1  N 

It, I  (6) 


It  should  be  noted  that  all  of  the  roughness  statistics  are  calculated 


using  the  centerline  as  the  datum  for  y.  This  is  defined  as  the 
datum  at  which  the  average  value  of  y  is  zero.  Rq  is  the  root  mean 
square  height  given  as: 


(V) 


R t  is  the  height  from  the  called  maximum  peak  to  the  minimum 
trough  and  is  given  as: 

Rt  Tmax  Tmin  (8) 

Rz  is  called  the  ten  point  height  and  is  given  as  the  mean  of  the 
difference  of  the  five  highest  peaks  and  the  five  lowest  troughs. 

_1  V 

^ z  T  (kmax,  Tmini)  (9) 

J  ,=  1 


The  correlation  length  scale,  \corr,  is  calculated  as  the  distance  (j 
times  the  sampling  interval,  Lp)  at  which  the  autocorrelation  func¬ 
tion  falls  below  0.5.  The  autocorrelation  function  is  given  by: 


N-  1- 


Cj=~ 


j  N-  1  -j 

1 — r  2  (y&i+j) 

1  j  i=  i 


i 

— y 

N-1& 


(10) 


It  should  be  noted  that  for  the  120-grit  and  smoother  surfaces  this 
value  was  less  than  the  sampling  interval  so  no  accurate  estimate 
could  be  made.  The  root  mean  square  slope,  slrms ,  is  given  as 
follows: 


si 


rms 


i  -y,  (Ti+t  yd 

FT M  u, ■+,-*,) 


2 


(ID 


A  similar  parameter,  the  root  mean  square  slope  angle,  was  of¬ 
fered  by  Acharya  et  al.  [12]  as  an  important  one  in  describing 
roughness  caused  by  machining  on  turbine  blades.  By  calculating 
the  power  spectral  density  of  the  surface  waveforms  using  a  fast 
Fourier  transform  and  the  first  three  even  moments  of  the  power 
spectral  density,  Townsin’s  height  parameter,  h,  was  calculated 
using  Eq.  (4).  Musker  [8]  offers  an  alternative  roughness  length 
scale  given  by: 

h'=Rq(l+aSp)(l+bSkKu)  (12) 

His  results  show  that  this  roughness  length  scale  works  well  for 
correlating  the  roughness  function  for  a  range  of  ship  hull  coatings 
when  a  long  wavelength  cutoff  of  2  mm  is  used  and  the  constants 
a  and  b  are  taken  to  be  0.5  and  0.2,  respectively. 


Uncertainty  Estimates 

Precision  uncertainty  estimates  for  the  resistance  measurements 
were  made  through  repeatability  tests  using  the  procedure  given 
by  Moffat  [23],  Ten  replicate  experiments  were  made  with  each  of 
the  test  plates  at  each  Reynolds  number.  This  was  carried  out  so 
that  the  relatively  small  differences  in  the  frictional  resistance 
between  the  surface  conditions  could  be  identified.  The  standard 
error  for  CF  was  then  calculated.  In  order  to  estimate  the  95% 
precision  confidence  limits  for  a  mean  statistic,  the  standard  error 
was  multiplied  by  the  two-tailed  t  value  (/=  2.262)  for  9  degrees 
of  freedom,  as  given  by  Coleman  and  Steele  [24],  The  resulting 
precision  uncertainties  in  CF  were  =S±0.3%  for  all  the  tests.  The 
overall  precision  and  bias  error  was  dominated  by  the  systematic 
error  due  to  the  combined  bias  of  the  force  gages  (±0.25%  full 
scale).  The  resulting  overall  precision  and  bias  uncertainty  in  CF 
ranged  from  ±4.8%  at  the  lowest  Reynolds  number  to  ±1.4%  at 
the  highest  Reynolds  number.  Periodically  throughout  the  experi¬ 
ments,  a  reference  plate  was  run  to  check  that  the  resulting  mean 
CF  value  was  within  the  precision  uncertainty  bounds  that  had 
been  obtained  from  previous  testing  with  the  same  surface.  This 
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Table  1  Roughness  statistics 


Specimen 

Af£, 

Ra 

R, 

R, 

Rz 

sk 

Ku 

Aeon- 

Slrms 

h 

(mm) 

Oun) 

(Mm) 

(Mm) 

(Mm) 

(Mm) 

(Mm) 

NF 

2.7 

5.0 

39 

38 

2.0 

5.0 

127 

0.113 

2.2 

±0.2 

±0.3 

±2 

±2 

±0.2 

±0.3 

±22 

±0.006 

+  0.3 

Unsanded 

10 

3.4 

4.8 

50 

43 

0.02 

4.8 

118 

0.113 

2.1 

±0.2 

±0.3 

±3 

±3 

±0.1 

±0.3 

±9 

±0.006 

±0.3 

2.5 

3.4 

28 

25 

-0.01 

3.4 

31 

0.112 

1.0 

±0.1 

±0.2 

±2 

±2 

±0.1 

±0.2 

±  1.7 

±0.006 

±0.1 

NF 

0.96 

1.63 

12 

II 

1.1 

3.2 

40 

0.063 

0.24 

±0.06 

±0.09 

±1.6 

±0.9 

±0.1 

±0.4 

±7 

±0.002 

±0.01 

60-Grit 

10 

1.14 

1.62 

13 

II 

0.61 

2.1 

31 

0.063 

0.24 

±0.07 

±0.08 

±  1.6 

±0.9 

±0.1 

±0.5 

±6 

±0.002 

±0.01 

1.05 

1.43 

13 

11 

-0.21 

1.5 

<25* 

0.063 

0.19 

1 

±0.05 

±0.07 

±2 

±0.6 

±0.1 

±0.2 

±0.002 

±0.01 

NF 

0.58 

1.02 

9.4 

7.6 

1.4 

5.2 

<25* 

0.047 

0.10 

±0.03 

±0.04 

+  0.7 

±0.2 

±0.2 

±  l.l 

±0.002 

±0.005 

120-Grit 

10 

0.70 

±0.02 

1.01 

±0.03 

10 

±0.95 

8.1 

±0.3 

0.96 

±0.09 

4.0 
±  1.0 

<25* 

0.047 

±0.002 

0.10 

±0.005 

0.65 

0.90 

9.2 

7.7 

0.11 

2.5 

<25* 

0.047 

0.08 

±0.03 

±0.03 

±0.6 

±0.2 

±0.08 

±0.6 

±0.002 

±0.005 

NF 

0.47 

0.86 

8.6 

7.0 

1.8 

7.8 

<25* 

0.039 

0.07 

±0.01 

+  0.04 

±  1.5 

±0.7 

±0.2 

±  1.9 

±0.002 

±0.006 

220-Grit 

10 

0.59 

±0.02 

0.86 

±0.04 

9.4 
±  1.6 

7.7 

±0.8 

1.1 

±0.2 

5.3 
±  1.3 

<25* 

0.039 

±0.002 

0.07 

±0.006 

0.52 

0.74 

9.1 

7.0 

0.15 

3.9 

<25* 

0.039 

0.05 

±0.02 

±0.04 

±  1.2 

±0.7 

±0.06 

±0.8 

±0.002 

±0.005 

NF 

0.43 

0.77 

8.0 

7.0 

1.8 

8.1 

<25* 

0.038 

0,06 

±0.03 

±0.04 

±0.5 

±0.7 

±0.1 

±  1.6 

±0.001 

±0.007 

400-Grit 

10 

0.52 

±0.03 

0.77 

±0.04 

8.6 

±0.4 

7.7 

±0.8 

1.4 

±0.2 

6.5 
±  1.3 

<25* 

0.038 

±0.001 

0.06 

±0.007 

0.49 

0.70 

8.1 

7.0 

0.5 

4.3 

<25* 

0.038 

0.05 

1 

±0.03 

±0.03 

±0.4 

±0.7 

±0.2 

±0.9 

±0.001 

±0.003 

NF 

0.40 

0.73 

7.9 

6.0 

2.1 

15 

<25* 

0.038 

0.05 

±0.01 

±0.02 

±2.6 

±0.6 

±0.7 

±  14 

±0.002 

±0.003 

600-Grit 

10 

0.47 

0.73 

8.4 

6.3 

1.8 

13 

<25* 

0.038 

0.05 

±0.02 

±0.03 

±2.5 

±0.6 

±0.7 

±  12 

±0.002 

±0.003 

0.47 

0.67 

9.1 

6.5 

0.5 

7.7 

<25* 

0.038 

0.04 

±0.02 

±0.03 

±2.5 

±0.9 

±0.3 

±6 

±0.002 

±0.003 

NF 

0.18 

0.30 

2.0 

1.9 

1.15 

3.3 

<25* 

0.015 

0.009 

+  0.004 

±0.01 

±0.3 

±0.3 

±0.08 

±0.3 

±0.001 

±0.001 

Polished 

10 

0.21 

±0.006 

0.30 

±0.01 

2.4 

±0.1 

2.1 

±0.4 

0.70 

±0.07 

2.3 

±0.2 

<25* 

0.015 
±  0.001 

0.009 

±0.001 

0.19 

0.27 

2.6 

2.2 

0.18 

2.1 

<25* 

0.015 

0.007 

±0.004 

±0.01 

±0.1 

±0.1 

±0.06 

±0.2 

±0.001 

±0.001 

*  -  Kon  was  <  Z* 

All  uncertainties  represent  the  95%  confidence  precision  bounds  for  the  measurement 


was  confirmed  in  all  cases  tested.  Uncertainty  estimates  for  the 
roughness  statistics  were  calculated  in  the  same  manner  and  are 
reported  in  Table  1. 

Results  and  Discussion 

The  presentation  of  the  results  and  discussion  will  be  organized 
as  follows.  First,  a  qualitative  discussion  of  the  nature  of  each  of 
the  surfaces  tested  will  be  made.  The  roughness  statistics  will  then 
be  presented.  Next  the  results  of  the  hydrodynamic  tests  will  be 
presented  and  discussed.  Finally,  an  attempt  will  be  made  to  relate 
the  roughness  statistics  of  this  class  of  surfaces  to  the  roughness 
function,  A  U+ . 

Qualitative  Description  of  the  Surfaces.  In  order  to  better 
understand  the  nature  of  each  of  the  surfaces  tested,  a  qualitative 
description  of  each  will  be  made  using  the  three-dimensional  to¬ 
pographic  profiles  shown  in  Fig.  2.  Even  cursory  inspection  of  the 
profiles  shows  that  the  surfaces  vary  greatly.  Figure  2(a)  shows 
the  unsanded  surface  and  indicates  that  it  has  relatively  large  fea¬ 
tures  with  a  wavelength  of  up  to  1  mm.  This  is  very  common  in 
as-sprayed  paint  surfaces  and  is  often  referred  to  as  “orange  peel” 
because  of  the  characteristic  texture.  Figure  2(b) ,  which  shows  the 
60-grit  surface,  indicates  that  the  orange  peel  has  been  almost 
entirely  removed  by  sanding,  but  linear  scratches  have  been 
added.  These  scratches  have  a  width  of  up  to  150  gm  and  a  depth 
of  up  to  25  /un.  Figure  2(c)  shows  the  120-grit  surface.  Many  of 
the  scratches  seen  in  the  60-grit  surface  have  been  removed,  and 
narrower,  shallower  scratches  have  been  added.  It  is  of  note  that 
some  of  the  deeper  scratches  from  the  60-grit  surface  have  not 
been  completely  removed  and  remain  as  features  of  up  to  150  /mi 
in  width,  with  a  depth  of  up  to  10  /mi.  The  220-grit  surface  (Fig. 
2(d))  shows  that  the  scratches  from  the  60-grit  paper  have  been 
removed.  They  have  been  replaced  with  finer  scale  scratches  that 
are  much  narrower  and  shallower.  The  400-grit  and  600-grit  sur¬ 
faces  are  very  similar  in  nature  and  for  this  reason  only  the  600- 
grit  surface  (Fig.  2(e))  is  presented.  The  polished  surface  (Fig. 
2(f))  shows  that  many  of  the  small  scale  peaks  and  troughs  seen 
in  the  400-grit  and  600-grit  surfaces  have  been  removed. 


The  evolution  of  the  surface  scratches  can  be  more  easily  seen 
in  plan  views  of  the  previous  figures.  These  views  for  the  un¬ 
sanded,  60-grit,  120-grit,  and  220-grit  surfaces  are  given  in  Fig.  3. 
Figure  3(a)  of  the  unsanded  specimen  shows  the  orange  peel  sur¬ 
face.  The  60-grit  surface  (Fig.  3(b))  shows  that  the  orange  peel  has 
been  removed  and  linear  scratches  have  been  added.  Smaller  scale 
scratches  are  evident  in  the  120-grit  surface  (Fig.  3(c))  as  well  as 
60-grit  scratches  that  have  not  been  completely  removed.  By  the 
time  the  220-grit  paper  has  been  used  (Fig.  3(d)),  only  rather 
small  scale  features  remain. 

The  quantitative  statistics  of  the  roughness  surfaces  are  given  in 
Table  1.  The  results  are  presented  for  processing  with  no  long 
wavelength  filter,  a  10  mm  long  wavelength  filter,  and  a  1  mm 
long  wavelength  filter.  One  item  of  note  is  that  all  of  the  rough¬ 
ness  tested  in  the  present  study  is  quite  small  compared  to  the 
roughness  used  in  a  majority  of  previous  studies.  Most  basic  re¬ 
search  has  focused  on  roughness  large  enough  to  generate  turbu¬ 
lent  flows  in  the  fully-rough  regime,  and  the  studies  on  ship  hull 
roughness  by  Musker  [8]  and  others  addressed  smaller  scale  tran¬ 
sitional  roughness  with  150  /rm^Rf=s600  /im  (5=Sfe+^320;  k+ 
based  on  R,).  For  the  present  study,  the  range  of  roughness  was 
2  /ims£R,5£39  /im  (0.15=£k+=£5;  k+  based  on  Rt).  This  is  im¬ 
portant  to  keep  in  mind  because  the  differences  in  drag  for  the 
surfaces  is  expected  to  be  rather  small.  The  change  in  the  rough¬ 
ness  statistics  for  the  unfiltered  profiles  with  sanding  is  shown  in 
Fig.  4.  The  figure  shows  that  all  of  the  roughness  height  param¬ 
eters  are  reduced  with  sanding  up  to  220-grit.  At  that  point,  no 
significant  reduction  in  the  roughness  height  is  made  by  sanding 
with  400-grit  and  600-grit.  The  polished  surface  does  show  a  sig¬ 
nificant  reduction  in  roughness.  Figure  2(f)  shows  that  this  is 
largely  due  to  a  reduction  in  the  isolated  protuberances  seen  in  the 
220-grit,  400-grit,  and  600-grit  surfaces. 

Test  Results.  The  results  of  the  hydrodynamic  tests  are 
shown  in  Fig.  5.  The  Schoenherr  mean  line  for  smooth  plates  is 
shown  for  comparison  [25].  The  Schoenherr  mean  line  is  given  as: 

0.242 

-=  =  log(ReLCV)  (13) 

vcF 

Table  2  shows  the  %  increase  in  CF  for  the  test  surfaces  com¬ 
pared  to  the  polished  surface.  The  results  show  that  the  60-grit 
specimen  had  a  significant  reduction  in  CF  from  the  unsanded 
surface.  A  further  reduction  in  CF  was  found  for  the  120-grit 
specimen  compared  with  the  60-grit  surface.  A  smaller,  but  sig¬ 
nificant,  reduction  in  CF  occurred  for  the  220-grit  surface  com¬ 
pared  with  the  120-grit  specimen.  However,  no  significant  change 
in  CF  was  seen  for  the  400-grit  and  600-grit  surfaces  compared  to 
the  220-grit  surface.  Inspection  of  Figs.  2(d-2e)  and  Table  1 
shows  that  the  roughness  on  these  surfaces  is  quite  similar  as  well. 
The  polished  surface  had  CF  values  that  were  significantly  less 
than  any  of  the  other  specimens.  The  reduction  in  frictional  resis¬ 
tance  seems  to  be  due  to  a  reduction  in  the  isolated  protuberances 
(Fig.  2(f))  that  were  seen  in  the  600-grit  surface  (Fig.  2(e)). 

It  should  be  noted  that  the  results  could  have  been  affected  to 
varying  degree  by  the  influence  of  the  surface  roughness  and  Rey¬ 
nolds  number  on  the  transition  of  the  flow  to  turbulence.  The  flow 
was  not  tripped,  so  the  transition  point  may  have  varied  between 
the  test  cases  depending  on  the  surface  roughness.  Leading-edge 
and  three-dimensional  effects  could  have  also  had  some  influence 
on  these  results.  However,  these  effects  are  very  difficult  to  quan¬ 
tify  precisely.  Overall,  it  is  felt  these  effects  are  small  since  the 
smooth  plate  results  agree  within  ~1%  with  the  Schoenherr  mean 
line  for  frictional  resistance  in  turbulent  flow.  The  fact  that  the 
smooth  plate  results  remain  a  fairly  constant  percentage  lower 
than  the  Schoenherr  mean  line  over  the  range  of  Reynolds  number 
tested  seems  to  also  indicate  that  the  transition  point  must  not  vary 
significantly  over  the  Reynolds  number  range  tested. 
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Fig.  2  Surface  waveforms  for  (a)  the  unsanded  specimen,  ( b )  the  60-grit  specimen,  (c)  the  120-grit  specimen,  (of)  the  220-grit 
specimen,  (e)  the  600-grit  specimen,  and  ( f )  the  polished  specimen.  (Uncertainty  in  the  y-direction  ±1  /um,  x-  and  z-directions  ±5 
/*m) 
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Fig.  3  Plan  view  of  the  surface  waveform  for  (a)  the  unsanded  specimen,  (b)  the  60-grit  specimen,  (c)  the  120-grit  specimen, 
and  ( d)  the  220-grit  specimen.  (Uncertainty  in  the  y-direction  ±1  /an,  x-  and  z-directions  ±5  /an) 


Specimen 


Fig.  4  The  effect  of  sanding  on  the  roughness  statistics  of  the 
unfiltered  profiles.  (Error  bars  represent  the  95%  confidence 
limits  for  the  precision  uncertainty.) 


Fig.  5  Overall  frictional  resistance  coefficient  versus  Rey¬ 
nolds  number  for  all  the  specimens.  (Precision  uncertainty 
^±0.3%  at  all  Reynolds  numbers;  overall  precision  and  bias 
ranges  from  ±1.4%  at  highest  Reynolds  number  to  ±4.8%  at 
lowest  Reynolds  number.) 
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Table  2  Increase  in  overall  frictional  resistance  coefficient  for 
the  test  specimens  compared  to  the  polished  surface 


Specimen 

Average  % 
Increase  in 

CF 

Range  of  % 
Increase  in 

CF 

Unsanded 

5.0 

3.0 -7.3 

60-Grit 

2.6 

1 

O 

c4 

120-Grit 

1.9 

1.0 -3.4 

220-Grit 

1.2 

0.8  -  1.7 

400-Grit 

1.2 

o 

i 

Uj 

600-Grit 

1.0 

i 

SO 

O 

Relation  of  Roughness  Statistics  to  \U+.  Using  similarity 
law  analysis,  Granville  [16]  derived  an  expression  to  relate  the 
local  frictional  coefficient,  Cf,  at  the  trailing  edge  of  a  planar 
surface  to  the  overall  frictional  resistance  coefficient,  CF,  for  the 
same  surface.  It  is  given  as: 


By  solving  this  equation  for  UT,  the  viscous  length  scale,  vlUT, 
at  the  trailing  edge  of  the  plate  can  be  obtained.  For  the  present 
surfaces  vlUT  ranged  from  ~14  /r.m  at  ReL=2.8X10s  to  ~7.6 
jum  at  ReL=5.5X  106.  The  roughness  function,  A U+,  at  the  trail¬ 
ing  edge  of  the  plate  can  be  found  using  the  method  of  Granville 
[16]  as  well.  This  procedure  involves  plotting  <J2/CF  versus 
ReLCF.  The  roughness  function,  A  U+ ,  is  given  as  the  following 
evaluated  at  the  same  value  of  ReL  CF  for  both  smooth  and  rough 
walls: 


In  the  present  study,  the  results  for  the  polished  surface  were 
used  as  smooth  plate  values.  Since  the  behavior  of  A  U+  at  van¬ 
ishing  roughness  height  did  not  behave  as  a  Nikuradse-type 
roughness  function,  attempts  were  made  to  collapse  the  results  to 
a  Colebrook-type  roughness  function  as  given  by: 

AI/+  =  — ln(l  +  k+)  (16) 

K 

All  of  the  roughness  length  scales  in  Table  1  were  considered 
including  the  Musker  [8]  and  Townsin  [20]  length  scales.  The  best 
fit  of  the  results  to  Eq.  (16)  was  obtained  using  a  multiple  of  the 
centerline  average  height,  Ru  ,  calculated  from  the  unfiltered  pro¬ 
files  as  k.  With  k—  1.35Ra  ,  87%  of  the  variance  (i.e.,  R2  =  0.87)  in 
A  U+  could  be  explained  with  the  Colebrook-type  roughness 
function  (Eq.  (16)).  The  results  are  shown  in  Fig.  6.  It  should  be 
noted  that  >80%  of  the  variance  could  also  be  explained  for  this 
relatively  simple  roughness  using  Rq  ,  R, ,  or  R,  calculated  from 
the  unfiltered  profiles.  Attempts  to  use  the  filtered  profile  statistics 
led  to  larger  scatter  in  the  roughness  function  than  the  unfiltered 
profile  statistics.  This  seems  to  indicate  that  for  sanded  surfaces  in 
this  Reynolds  number  range,  long  wavelength  roughness  (up  to  50 
mm)  contributes  significantly  to  the  increase  in  frictional  resis¬ 
tance. 

Using  the  roughness  function  obtained  for  a  flat  plate,  Granville 
[16]  gives  a  similarity  law  procedure  for  calculating  the  effect  of 


Fig.  6  Roughness  function  for  all  specimens.  (Overall  uncer¬ 
tainty  ±0.1  in  A  U+) 


the  same  roughness  on  a  planar  surface  of  arbitrary  length.  This 
was  carried  out  with  the  present  results  for  surfaces  of  4  m  and  12 
m.  These  lengths  were  chosen  to  be  representative  of  the  range  of 
length  of  sailing  vessels  for  which  sanding  would  be  practicable. 
Using  this  analysis,  significant  increases  in  CF  above  the  polished 
plate  values  were  predicted  for  the  surfaces.  These  increases  were 
of  the  same  order  as  seen  in  the  present  experiments.  For  example, 
it  was  found  that  the  unsanded  surface  would  have  an  average 
increase  in  CF  of  5.0%  above  a  polished  surface  over  a  velocity 
range  of  2. 3-4. 6  m/s  (4.5-8. 9  knots)  at  a  length  of  4  m  and  an 
average  increase  in  CF  of  4.5%  over  a  velocity  range  of  2. 6-5.0 
m/s  (5. 1-9.7  knots)  at  a  length  of  12  m.  Over  the  same  velocity 
range,  the  60-grit  surface  would  have  an  average  increase  in  CF  of 
2.5%  at  a  length  of  4  m  and  an  average  increase  in  CF  of  2.3%  at 
a  length  of  12  m.  The  increase  for  the  120-grit  surface  would  be 
1.9%  at  a  length  of  4  m  and  1.6%  at  a  length  of  12  m.  Both  the 
220-grit  and  400-grit  surfaces  would  have  an  increase  of  1 .3%  and 
1.1%  at  lengths  of  4  m  and  12  m,  respectively,  while  the  increase 
for  the  600-grit  surface  would  be  1.1%  and  1.0%  at  lengths  of  4  m 
and  12  m,  respectively.  Due  to  the  three-dimensional  nature  of  the 
boundary  layer  flow  around  a  real  vessel,  however,  additional  data 
on  actual  hull  shapes  are  needed  to  corroborate  these  findings.  It 
should  also  be  noted  that  in  the  present  study  significant  care  was 
given  to  ensure  that  the  entire  surface  was  sanded  and  that  no  area 
was  missed.  This  was  relatively  easy  on  a  small  flat  surface,  but 
would  be  much  more  difficult  in  practice  on  a  large  three- 
dimensional  shape. 


Conclusion 

Measurements  of  the  roughness  and  frictional  resistance  of 
sanded  paint  surfaces  have  been  made.  The  results  indicate  that 
as-sprayed,  unsanded  surfaces  can  have  a  significant  increase  in 
CF  compared  to  a  polished  surface.  Smaller,  but  significant,  in¬ 
creases  in  CF  compared  to  the  polished  surface  were  also  noted  on 
surfaces  sanded  with  sandpaper  as  fine  as  600-grit.  This  increase 
seems  to  be  due  to  isolated  surface  protuberances  not  completely 
removed  by  the  sanding  process.  The  roughness  function,  AI/+, 
shows  good  collapse  to  a  Colebrook-type  roughness  function  for 
this  class  of  surfaces  when  a  multiple  of  the  centerline  average 
height  (fe=1.35R„)  is  used  as  the  roughness  length  scale.  Simi¬ 
larity  law  predictions  of  CF  on  larger  planar  surfaces  of  sailing 
vessel  length  show  that  similar  increases  in  CF  can  be  expected  in 
that  range  of  wetted  length  as  well.  Further  effort  needs  to  be 
focused  on  understanding  the  effect  of  roughness  on  three- 
dimensional  bodies. 
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Nomenclature 

a,  b  =  constants  in  Musker’s  roughness  length  scale  equa¬ 
tion 

B  =  smooth  wall  log-law  intercept=5.0 
CF  —  overall  frictional  resistance 
coefficient  =  (FD)/(  1/2  pU2eS) 

Cf  =  local  frictional  resistance  coefficient=(r0)/(l/2pUp) 

C j  =  autocorrelation  function 
E  =  power  spectral  density  of  surface  waveform 
Fd  —  drag  force 

h  —  Townsin’s  roughness  height  parameter=  yjam0m1 
j  —  lag  in  autocorrelation  function 
k  =  arbitrary  measure  of  roughness  height 
ks  =  sand  roughness  height  or  equivalent  sand  roughness 
height 

Ku  —  kurtosis 
L  =  plate  length 

Lp  =  overall  length  of  surface  profile 
Ls  =  profile  sampling  interval 
[/]  =  other  roughness  length  scales 
m„  =  nth  moment  of  the  power  spectral  density 
N  =  number  of  samples  in  surface  profile 
Re ,5*  =  displacement  thickness  Reynolds  number=  UeS*/v 
ReL  =  Reynolds  number  based  on  plate  length=  UeLlv 
Ra  —  centerline  average  roughness  height=  (1//V)2^1 j | y ,■  | 

Rq  —  root  mean  square  roughness  height=  \/(l //V)2^L  iyj 
R ,  =  maximum  peak  to  through  height =ymax—ymjn 
Rz  =  ten  point  roughness  height=2j=1(ymax!—ymini) 

S  —  wetted  surface  area 

slrms  ~  root  mean  square  slope  of  the  roughness  profile 

=  V(l//V-l)2f=-11{(yi+1-yi)/(xi+1-xi)}2 

Sk  —  skewness 

Sp  —  average  absolute  slope  of  roughness  profile 
U  =  mean  velocity  in  the  x  direction 
Ue  =  freestream  velocity  relative  to  surface 
A  U+  =  roughness  function 
U T  =  friction  velocity  =  Jr0/p 
x  =  streamwise  distance  from  plate  leading  edge 
y  =  normal  distance  from  the  boundary  measured  from 
roughness  centerline 
=  bandwidth  parameter=n70m4/m2 
8  =  boundary  layer  thickness 
8*  =  displacement  thickness  =  Jq(1  —  U/Ue)dy 
k  —  von  Karman  constant  =0.41 
y  =  wavenumber=2i7-/\ 

X  =  wavelength 
XCOrr  =  correlation  length  scale 

v  —  kinematic  viscosity  of  the  fluid 


II  =  wake  parameter 
p  =  density  of  the  fluid 
ra  =  wall  shear  stress 
a)  =  wake  function 

Superscript 

+  =  inner  variable  (normalized  with  UT  or  UT/v) 

Subscript 

FL  =  long  wavelength  filter 
min  =  minimum  value 
max  =  maximum  value 
R  =  rough  surface 
S  =  smooth  surface 
TE  =  trailing  edge 
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